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ABSTRACT 

The age of the Universe has been increasingly constrained by different techniques such as the 
observations of type la supernovae at high redshift or dating the stellar populations of globular 
clusters. In this paper we present a complementary approach using the colours of the brightest 
elliptical galaxies in clusters over a wide redshift range (z ^ 1). We put new and independent 
bounds on the dark energy equation of state parameterised by a constant pressure-to-density 
ratio wq as well as by a parameter which determines the scaling between the matter 
and dark energy densities. We find that accurate estimates of the metallicities of the stellar 
populations in moderate and high redshift cluster galaxies can pose stringent constraints on 
the parameters which describe the dark energy equation of state. Our results are in good 
agreement with the analysis of dark energy models using SNIa data as constraint. Accurate 
estimates of the metallicities of stellar populations in cluster galaxies at z 2 will make this 
approach a powerful complement to studies of cosmological parameters using high redshift 
SNIa. 

Key words: galaxies: elliptical and lenticular, cD — galaxies: evolution — cosmology: Cos- 
mic Microwave Background, anisotropy, power spectrum 



1 INTRODUCTION 

The discovery that the evolution of the Universe may be dominated 
by an effective cosmological constant is one of the most remarkable 
findings of recent years. An exceptional opportunity is now opening 
up to decipher the nature of dark energy, to test the veracity of 
theories and to reconstruct the dark energy's profile using a wide 
variety of observations over a broad redshift range. 

One energy candidate that could possibly explain the obser- 
vations is a dynamical scalar "quintessence" field. In many cases 
these models alleviate the fine-tuning and coincidence problems 
carried by a cosmological constant (Zlatev et al. 1999). The most 
relevant observational evidence for quintessence is that its equa- 
tion of state parameter, wq ~ Pq/pq can assume values different 
from —1, which is the expected value in the case of a cosmolog- 
ical constant. Moreover, since very little is known currently about 
the nature of dark energy, an interaction (aside from gravity) of 
quintessence with the remaining forms of matter could be expected 
as well (e.g. Wetterich 1995). 

Other deviations from a cosmological constant can be con- 
sidered as, for example, a sound speed different from light or 
perturbations in the dark energy fluid. However these aspects of 
quintessence affect mainly the growth of perturbations and their 
detectability is therefore subjected to the assumption of a model 
of structure formation. In this work we will be most interested in 
the effect that dark energy causes on the background evolution. In 
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particular, we show that an accurate determination of the age of the 
Universe can be used to obtain important and independent informa- 
tion on the dark energy component. 

Discrepancies between age determinations have long plagued 
cosmology. The situation has changed drastically in recent years: 
Type la supernova measurements (Perlmutter et al. 1999; Riess et 
al. 1998), the acoustic peaks in the CMB anisotropics (Netterfield 
et al. 2001), cosmochronology of long-lived radioactive nuclei in 
old Population II stars (Cayrel et al. 2001), just to name a few, are 
all consistent with an age to = 14 ± 1 Gyr. Furthermore, the recent 
age determination of the oldest white dwarf in the globular cluster 
M4 by the Hubble Space Telescope, with U,d = 12.0 ± 0.7 at 95% 
c.l. (Hansen et al. 2002), can be used to put the strongest lower limit 
to date on the age of the Universe. 

It is therefore timely to investigate if the new age determina- 
tions can be used to produce new constraints on dark energy. These 
results could also prove very useful in breaking some of the degen- 
eracies that arise in cosmological parameter estimations from other 
observational sources (e.g. Bean & Melchiorri 2002; Maor et al. 
2002). 

Recently, Krauss & Chaboyer (2001), have provided con- 
straints on Wq by using new globular cluster age estimates. Our 
analysis will differ in two new fundamental aspects: first we will 
use age estimates based on old stellar populations at moderate red- 
shift. As in (Ferreras, Melchiorri & Silk 2001) we will provide 
estimates on the age of the Universe via a comparison of pop- 
ulation synthesis models with the photometry of the oldest stel- 
lar populations in bright cluster ellipticals. Second, following the 
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parametrization of Dalai et al. (2001) we will also investigate the 
possibility of an interaction besides gravity between dark energy 
and dark matter 



2 MODELS OF DARK ENERGY AND THE AGE OF THE 
UNIVERSE. 

In the literature one can find many proposals of scalar field dark en- 
ergy models, differing by the choices of the self-interacting poten- 
tials (e.g. Sahni 2002), resulting in equations of state that in general 
depend on time. On the other hand, the age of the Universe is most 
influenced by the relatively late evolution history. This enables us 
to simplify the problem by assuming a constant equation of state 
for the dark energy. A parameter determination through a likeli- 
hood analysis would very probably have a few chances of detecting 
a relatively smooth variation of the equation of state. This line of 
reasoning is very similar to the one usually expressed in many stud- 
ies of constraint from SN-Ia (Maor et al. 2002). 

In fact, as our results will show, the SNe parameter degener- 
acy is significantly similar to our case. Thus our work in a way 
offers an independent check of those findings. This has its own 
relevance due to the intense debate on whether all the important 
systematic errors have been taken into account in SNe data reports. 
In what follows the sub-indices Q, C and B indicate respectively 
the contributions from dark energy-quintessence, cold dark matter 
(CDM) and baryons. The collective CDM and baryonic contribute 
will be represented by the label Al. We will adopt a quite general 
phenomenological model in which the ratio of the dark energy and 
dark matter densities scale as a power law: 

PQ oc pca^, 



(1) 



with ^ kept constant (as in Dalai et al. 2001). It should be noted 
that this assumption generally involves an interaction besides grav- 
ity. The dark energy equation of state, wq = Pq/pq, is held fixed 
as stated above. Three interesting cases deserve particular atten- 
tion. Putting ^ = 3 and wq = — 1 gives the ACDM model. The 
condition ^ = —3wq corresponds to decoupled quintessence. Fur- 
thermore the choice ^ = implies models that may explain the 
coincidence problem through an identical scaling of the two energy 
densities (e.g. Tocchini-Valentini & Amendola 2002, Amendola & 
Tocchini-Valentini 2001). The baryons are assumed coupled only 
through gravity to the rest of the world. 

To perform a comparison with the data we will need an ex- 
pression of the age of the Universe for this class of models. This is 
given by: 



to = 9.8Gyr h' 
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Figure 1. Evolution with age of the colour of bright cluster galaxies at 
moderate redshift. The shading corresponds to metallicities in the range 
0.0 < < 0.4 dex. The dashed lines represent the observed colour range 
of the brightest galaxies in the three clusters presented. The arrows mark 
the allowed range at the observed redshift given the colour constraint (see 
text for details). 



The inclusion of baryons is important since their contribution may 
become significant not too far in the past history of the Universe 
and their exclusion may overestimate the age of the Universe by 
more than 10% relative to the complete case, especially for highly 
negative values of wq. 

In the decoupled case it has been shown that values of the 
equation of state such that wq < —1 are not completely ex- 
cluded by observations (Melchiorri et al 2002, Schuecker et al. 
2002). Therefore in our analysis we will consider the interval 
—3 < Wq < 0. The limitations < ^ < 9 will be imposed on 
the parameter ^. To reduce the number of free parameters we con- 
sider a flat Universe. For models of coupled quintessence we fix 
J^jvi ~ 0.3 for the density of CDM and baryons altogether, and we 
also fix Sis = 0.04 for baryons. 



3 USING BRIGHT ELLIPTICALS AS DARK ENERGY 
COSMIC CLOCKS 

The use of bright cluster ellipticals over a wide range of redshifts 
allow us to set contraints on the age of the Universe. In the volume 
sparmed by the various cosmological parameters involved, these 
constraints go roughly "in the same direction" as the constraints 
from high redshift SNe data — which depend on the luminosity 
distance, and thus provide an independent method to explore the 
allowed values of the parameters. Ferreras et al. (2001) used this 
approach to determine the age of the Universe and showed that the 
results were consistent with estimates using other observations such 
as high redshift SNIa (Riess et al. 1998; Perlmutter et al. 1999), age 
dating of globular clusters (Salaris & Weiss 1998; Pont et al. 1998; 
Carretta et al. 2000) or nucleocosmochronology of old halo stars 
(Cayrel et al. 2001). 

FigureQshows the methodology of this approach. The colours 
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Figure 2. Parameter space allowed by the colour constraints from the 
brightest ellipticals in the three clusters shown in figure □ Var lous cosmo- 
logical parameters are fixed, namely: Sljv/ = 0.3; Hq = 72km s~^Mpc^^ 
and Qb = 0.04. The allowed region of parameter space stays to the right 
of the shaded line. We have chosen two of the latest formation redshifts 
compatible with the observations: zp = 2 (top panel) and zp = 3. Higher 
values of 2^ will further constrain the allowed area of parameter space to 
the right. 



of the first few brightest ellipticals from clusters at moderate red- 
shift is compared with the predictions of the population synthe- 
sis models of Bruzual & Chariot (1993) for a range of ages and 
metallicities. The latest available version of these models have been 
used. The shaded area represents the colours for the metallicity 
range < < +0.4 dex, as expected in the central regions 
of bright ellipticals (Henry & Worthey 1999). The dashed lines in 
each panel correspond to the observed colours of the brightest el- 
lipticals in clusters: CL 1358-1-62 (z = 0.33; Van Dokkum et al. 
1998); MS 1054-^03 {z = 0.83; Van Dokkum et al. 2000), and 
J0848-I-4453 (z = 1.27; Van Dokkum & Stanford 2002). The ar- 
rows give the allowed age of the Universe at the redshift of the clus- 
ter for the metallicity range expected in bright early- type galaxies. 
The upper limit can be used to constrain the allowed values of the 
parameters wq and ^. Figure |2| shows the constraint imposed by 
the upper limit expected for the three clusters shown in the previ- 
ous figure. Some of the cosmological parameters have been fixed, 
namely: fijv/ — 0.3, JIq = 1 — JIm, Ho ~ 72 km s~^ Mpc^^, 
and = 0.04. The allowed region stays to the right of the shaded 
line for each of the three clusters. As reference, a ACDM cos- 
mology corresponds to (wq,^) = ( — 1,3). A formation redshift 
of «F ~ 2 is assumed, an arguably late epoch for the formation 
of massive elliptical galaxies given the redshift evolution of the 
colour-magnitude relation (e.g. Stanford, Eisenhardt & Dickinson 
1998). Earlier formation epochs are possible and will further con- 
strain the region of parameter space to the right of the figure. 

A more comprehensive analysis was performed using the 
complete sample of 17 clusters from Stanford et al. (1998) over 
the redshift range 0.3 < z < 0.9. This sample was extracted 
from a larger sample of 46 clusters drawn from a variety of optical. 



X-ray and radio-selected clusters. The morphological classification 
was done from high resolution images taken by the Wide Field and 
Planetary Camera 2 on board the Hubble Space Telescope. We se- 
lected the three brightest early-type galaxies for each cluster, and 
used the published photometry in NIR passbands J, H and K, as 
well as on two optical filters which straddle the 4000A break. Along 
similar lines as the method shown above, we compare the photom- 
etry with the predictions of the latest population synthesis models 
of Bruzual & Chariot (1993). Furthermore, the stellar populations 
in these bright galaxies are assumed to evolve into those observed 
in bright ellipticals in local clusters such as Coma or Virgo. Hence, 
a further constraint is imposed comparing the U — V colour when 
the galaxy is evolved to zero redshift. 

The models are run assuming a flat cosmology. The dark en- 
ergy parameters are explored over a large range: —3 < wq < 0; 
and either < ^ < 9 or 0.1 < Qai < 1 for coupled and 
uncoupled models, respectively. The former assume fixed values 
for r2i\/ = 0.3 and Qb = 0.04. Uncoupled models are obvi- 
ously independent of the baryon density. We assume a wide range 
of formation redshifts: 2 < zp < 10 and stellar metallicities: 
-0.1 < + 0.3. A Salpeter (1955) initial mass function (IMF) 
was used, although the results - which are purely constrained by 
galaxy colour, not absolute luminosity - are rather insensitive to 
the choice of IMF. 



4 DISCUSSION 

The results of our analysis are reported in figures|3|and|4] For each 
theoretical model we compare the ages obtained through equation|2| 
with the observational data. Likelihood contours at the 68% (solid), 
and 95% (dashed) confidence levels are shown in the Q, m ~wq (un- 
coupled, Fig.|3} and ^ — wq (coupled, Fig.|4j planes. Since the re- 
sults depend also on 2 extra parameters involving the stellar popu- 
lations, namely the metallicity (Zt ) and the formation redshift (zp) 
we have investigated different possible values for these two param- 
eters. In all the plots we have assumed a gaussian distributed prior 
on the Hubble parameter ho — 0.72±0.08 as indicated by the Hub- 
ble Space Telescope Key project (Freedman et al. 2001). In Fig.|3| 
we report the results of our analysis assuming uncoupled dark en- 
ergy models and for different choices of Zi, and Zi?. As we can see 
from the plots a degeneracy is present in the tradeoff between wq 
and r2i\/ (higher values of \wq \ are consistent with a higher Q,m) 
and no strong independent constraint can be obtained on these pa- 
rameters. Furthermore, there is a strong and similar dependence on 
formation redshift and metallicity. Namely, an increase in zf or 
Zi, shifts the likelihood contours towards higher values of Qai. A 
model with r2j\/ = 1 and wq = 0, for example, is excluded at high 
significance when "reasonable" values of zp — 2 and — Zq 
are assumed, while an increase in formation redshift or metallic- 
ity pushes this model back to the 68% confidence level. On the 
other hand, a low value of can be used as a way to revert the 
argument: if we assume a "concordance" model with ^Im = 0.3 
and Wq = — 1, then high values of the formation redshift and the 
metallicity of the stellar populations in bright cluster ellipticals are 
disfavoured. Marginalizing over zp and Z^ gives the upper limits 
WQ < —0.35 and flu < 0.66 at 68% C.L.. These limits even 
if less stringent are in good agreement with those obtained from 
observations of luminosity distances of type la supernovae. 

In figure |4] we plot the contours in the ^ — wq planes fixing 
Qm = 0.3. The thick diagonal line represents uncoupled models. 
The figure shows that even if Qm is kept fixed, no strong limit on 
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Figure 3. Likelihood contours at tlie 68% (solid) and 95% (dashed) confi- 
dence levels in the parameter space spanned by Q,m and wq. Uncoupled 
dark energy models are considered. and zf represent the metallicity 
(with respect to solar) and the formation redshift of the stellar populations, 
respectively. 
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Figure 4. Likelihood contours at the 68% (solid) and 95% (dashed) confi- 
dence levels in the parameter space spanned by ^ and wq for coupled dark 
energy models. All models in this figure assume a fixed Q,m = 0.3. 
and zf represent the metallicity (with respect to solar) and the formation 
redshift of the stellar populations. Uncoupled models (g = — Sioq) are 
represented by the thick line. 



Wq can be placed. We cannot disentangle the degeneracy between 
WQ and ^ if low metallicities (Zi, ~ Zq) or late formation epochs 
{zf ~ 2) are assumed. However, coupled models 7^ — Siug) are 
not favoured in our analysis, unless one is willing to have high zf 
and Zi,. We should emphasize that 2f is treated as a free parameter 
that determines the ages of the stellar populations. Structure for- 
mation studies for the models explored in this paper would provide 
further constraints on the allowed region of parameter space. Nev- 
ertheless, the models with Zi, = IAZq are disfavoured in the light 
of the results from high redshift SNIa (Dalai et al. 2001) which rule 
out models with wq > —0.5 regardless of ^ at the 99% C.L. Even 
models with ^* = 1.2Zq would be hard to reconcile with the SNIa 
data if the formation redshift is too high. 

Figures |3| and |4] illustrate the fact that an accurate estimate 
of the metallicity along with a combined analysis with high red- 
shift SNIa would pose strong constraints on dark energy param- 
eters. High metallicities in massive ellipticals are expected if we 
assume a scenario with a highly efficient star formation in which 
most of the ejecta from stars is recycled back into subsequent gen- 



erations of stars. Wind models for elliptical galaxy formation sup- 
port a mass-metallicity relation which is caused by the potential 
well of the galaxy which prevents significant gas outflows in mas- 
sive galaxies thereby achieving high metallicities (Larson 1974; but 
see Silk 2003). Massive elliptical galaxies such as those targeted in 
this study should have solar or even supersolar average metallici- 
ties. However, the current data does not allow us to set a stringent 
constraint on Z*. Hence, in addition to advancing our knowledge 
of galaxy formation and evolution, accurate measurements of stel- 
lar metallicities in galaxies at moderate and high redshift will have 
a very significant impact on estimates of the cosmological parame- 
ters. 

The likelihood regions shown in the figure are in good agree- 
ment and/or consistent with those obtained by the analysis of high 
redshift SNIa (Dalai et al. 2001) for most of the models explored in 
this paper. Furthermore, our decoupled models (figure|3} agree well 
with the estimates from globular cluster ages (Krauss & Chaboyer 
2001). Hence, the study of the colours of the stellar populations in 
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moderate redshift elliptical galaxies deserves a place of its own as 
an estimator of cosmological parameters. 
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